Fermionization and fractional statistics 
in the strongly interacting one-dimensional Bose gas 
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We discuss recent results on the relation between the strongly interacting one-dimensional Bose 
gas and a gas of ideal particles obeying nonmutual generalized exclusion statistics (GES). The 
thermodynamic properties considered include the statistical profiles, the specific heat and local pair 
correlations. In the strong coupling limit 7 — + 00, the Tonks-Girardeau gas, the equivalence is 
with Fermi statistics. The deviation from Fermi statistics during boson fermionization for finite but 
large interaction strength 7 is described by the relation a ~ 1 — 2/7, where a is a measure of the 
GES. This gives a quantitative description of the fermionization process. In this sense the recent 
experimental measurement of local pair correlations in a ID Bose gas of 8r Rb atoms also provides a 
measure of the deviation of the GES parameter a away from the pure Fermi statistics value a = 1. 
Other thermodynamic properties, such as the distribution profiles and the specific heat, are also 
sensitive to the statistics. They also thus provide a way of exploring fractional statistics in the 
strongly interacting ID Bose gas. 

PACS numbers: 03.75.-b, 05.30.Jp, 05.30.Pr 



The one-dimensional (ID) Bose gas with delta-function 
interaction 0,0 is an elegantly simple exactly solved 
quantum many-body system which is testable in exper- 
iments on trapped quantum gases of ultracold atoms 
H II H H II The model has been extensively 
studied (see reviews [I3> EH E3, E3, E3> EM E3l ) The 
novel quantum many-body effects inherent in the ID in- 
teracting Bose gas continue to attract theoretical and ex- 
perimental interest [E 0, EE E3, 123, E3, EE HI, ISE E3l . 

In the experiments an effectively ID system of bosons 
is made by tightly confining the atomic gas in two radial 
directions and weakly confining it along the axial direc- 
tion. In this way the motion of particles along the radial 
directions is frozen out, thus making the system effec- 
tively ID. The delta-function interaction between bosonic 
atoms can be realized via short range interaction with an 
effective coupling constant g±u determined through an 
effective ID scattering length aiD [27j ■ In terms of the 
atomic mass m the dimcnsionlcss coupling constant 



7 



mgip 
h 2 n 



(1) 



is varied to characterize different interaction regimes by 
tuning the coupling strength giD or the particle number 
density n. In the weak coupling regime 7 <C 1, the ID 
Bose gas can undergo a quasi-Bose-Einstein condensation 
due to the coherence of individual particle wave func- 
tions. In the strong coupling regime 7 > 1, the bosons 
behave like impenetrable hard core particles, known as 
the Tonks-Girardeau gas (TG) 

The experimental realization of the ID TG gas [E IE IE 
0,0,0 has provided significant insight into the "fermion- 
ization" of interacting bosons. During the fermionization 
process the local pair correlations of strongly interacting 
bosons gradually vanish as the interaction strength giD 
increases [2E HE El HE H3- This behaviour was ob- 
served clearly in the local pair correlation function ob- 



tained from measurements of photoassociation rates in 
a ID Bose gas of 87 Rb atoms In particular, the 
wave functions of two colliding particles tolerate a de- 
gree of overlap for finitely strong interaction. This sug- 
gests that the behaviour of strongly interacting bosons 
deviates from pure Fermi behaviour for finite interaction 
strength. In this sense the quasiparticle excitations of 
strongly interacting bosons confined in ID do not obey 
pure Fermi statistics. 

In ID the dynamical interaction and the statistical in- 
teraction are inextricably related due to the pairwise in- 
teraction between particles. As a result ID interacting 
systems can map to a system of ideal particles obeying 
Haldane exclusion statistics [34Ll35| . Haldane formulated 
a general description of quantum statistics based on gen- 
eralized Pauli exclusion statistics, now called generalized 
exclusion statistics (GES). For GES the dimensionality of 
the single particle Hilbert space depends linearly on the 
particle numbers of other species when certain species of 
particles are added. Soon after the formulation of GES, 
Wu [3^1 found that the most probable distribution is de- 
termined by GES through counting the dimensionality 
of the single particle Hilbert space. This definition al- 
lows different species to refer to identical particles with 
different quasimomenta in a wide range of ID integrable 
models [35j. Over a decade ago it was shown that the ID 
Calogero-Sutherland model [3(| and the ID interacting 
Bose gas [IE H3 are equivalent to an ideal gas obeying 
GES. This equivalence has been recently investigated for 
a ID model of interacting anyons |38j . An anyon-fermion 
mapping has been proposed to obtain solutions for sev- 
eral models of ultracold gases with ID anyonic exchange 
symmetry [H). 

Here we discuss the statistical profiles and the thermo- 
dynamic properties of the strongly interacting ID Bose 
gas through the GES and thermodynamic Bethe ansatz 
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(TBA) jljj approaches. The results presented here are a 
special case of results [3^, |4l| obtained recently for the 
more general ID model of interacting anyons [381 li^ . 
which contains the ID interacting Bose gas as a spe- 
cial case. Here we highlight and discuss these results 
in the context of the strongly interacting ID Bose gas. 
At low temperatures the quantum statistics of strongly 
interacting bosons confined in ID is nonmutual, i.e., the 
GES parameter introduced in the next Section is given 
by c*ij = a 5jj , which also defines <?-on or fractional 
statistics |43|, |44|, |45| . The strongly interacting ID Bose 
gas implements a continuous range of GES, approaching 
Fermi statistics with a — 1 as the interaction strength 
gin — > ao. We discuss the total energy, specific heat and 
local pair correlations in terms of the interaction strength 
and the temperature. The results indicate that experi- 
ments on the strongly interacting ID Bose gas may be 
used as a testing ground for observing fractional statis- 
tics. 



I. BETHE ANSATZ AND GES 



The ID model of TV interacting bosons is defined by 
the Hamiltonian 0,0 
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under periodic boundary conditions. The coupling con- 
stant is determined by giD = h 2 c/m = 2h 2 / (m\aiu\) 
where the coupling strength c is tuned through an effec- 
tive ID scattering length oid via confinement Here- 
after we set H = 2m = 1 for convenience. We shall restore 
physical units in the thermodynamics later. The energy 
eigenvalues are given by E — Y^jLi kj> where the indi- 
vidual quasimomenta kj satisfy the Bethe equations [l| 
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for j = 1, . . . , N. Here L is the length of the system with 
n = N/L the particle density. 

In the thermodynamic limit, L —* oo with n finite, 
the ground state energy per particle is given by E/N — 
n 2 e(7) where 



p{x)x 2 dx. 



(4) 



The quasimomentum distribution function p(x) and the 
parameter A = c/Q, where Q is the cut-off momentum at 
zero temperature, are determined by the integral equa- 
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In the strong and weak coupling limits, the leading terms 
in the ground state energy can be obtained analytically 
from the Bethe equations © [H IHIH or © [1 l47fl . 
In the strong coupling regime 7 > 1 the ground state 
energy per particle is given by 
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The low-lying excitations close to the Fermi surface have 
a linear dispersion E = Eq + pv c with sound velocity 
v c = — 47 _1 ) as p — > in the thermodynamic limit. 
Here the Fermi velocity wp = 27m. Further, from the 
discrete Bethe equations (0 the finite-size corrections to 
the ground state energy for st rong coupling in the ther- 
modynamic limit are given by 46] 
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with central charge C = 1. 

Now consider the ID interacting Bose gas from the 
perspective of GES. The GES parameter a™ is defined 



through the linear relation Adi/ANj 



34], 



the number of available single particle states of species i, 
denoted by di, depends on the number of other species 
{Nj} when one particle of species i is added. Thus di is 



given by [35l 13 



d i ({N j })=G?-'£a ij N j 
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Here G° = rf({0}) is the number of available single par- 
ticle states with no particles present in the system. The 
GES parameter can be determined by the two-body scat- 
tering matrix 

oiij := a[k, k 1 ) = 5{k, k') - — 6'(k - k'), (9) 

27T 

through counting the available states in quasimomentum 
space. Here the function 



9'(x) = 



2c 



(10) 



All accessible states Di = di + Ni — 1 in a momentum in- 
terval Aki can be determined by the Bethe ansatz result 



P + Ph 
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dk'6'(k - k')p(k') 
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Here p and ph are the density of occupied states and the 
density of holes in interval Aki , respectively. 
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The equivalence between ID interacting bosons and 
ideal particles obeying GES is clearly seen from the 
equivalence between the TBA 0] result 



ik) = e°(k)-p 



T 
'2tt 



dk'9'(k - k') ln(l + C -< k '^ KBT ll2) 



and the GES equation |3J 

IV j 
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on setting 



e(ki)/K B T 



Here e(fc) is the dressed 



energy measuring the energy above the Fermi surface, 
with e°(k) = k 2 . 

In the above equations, the most probable distribution 
is given by J2j( w ^ig. + Pij) n j = 1 witn n i = N i/Gi 



and Pij = OijGj/Gi |35j. For a = the GES result O 



reduces to Bose statistics. For a; 



1 it reduces to Fermi 



statistics. The above equivalence is valid for the case 
c > 0. We note that this type of connection between TBA 
and GES is very general. For ID interacting systems with 
internal degrees of freedom the equivalence appears to be 
with ideal particles obeying multi-component exclusion 
statistics |ia El IHol Isil] . 



II. DISTRIBUTION PROFILES AND 
THERMODYNAMIC PROPERTIES 

In this section we concentrate on the distribution pro- 
files and the thermodynamics of the strongly interacting 
ID Bose gas. In general, the ID interacting many-body 
systems have mutual GES with quasiparticle excitations 
in momentum space. The quasiparticle excitations in 
these systems can be viewed as anyons with inperfect ex- 
clusion. In the GES description above, the state count- 
ing only involves real qusimomcntum roots. For complex 
quasimomentum roots, Wu's statistics (|13|) is not valid. 
For the ID Bose gas, we restrict our attention to the re- 
pulsive regime, i.e., gin > 0. From IT51) we see that at 
zero temperature there exists a Fermi-like surface with 
energy 6f = fM), where n(k) = l/a if e(k) < p,Q and 
n(k) = if e(k) > po- Therefore at low temperatures, a 
relatively small number of particles are excited above the 
Fermi surface leaving an unequal amount of holes below 
the Fermi surface. For strong coupling and at low tem- 
peratures, the majority of particles are below the surface. 
It is reasonable to expect that the statistics of the par- 
ticles at finite temperatures are equal to the statistics at 
zero temperature H3,0,0| . Varying the repulsive inter- 
action from weak coupling 7 = + to the strong coupling 
limit 7 — > 00 is equivalent to varying from Bose statistics 
a«0to Fermi statistics with a — > 1 for ideal particles 
obeying GES. We see from the TBA result JT^J) that for 
c = 0+ 



e(k) 



fj, 
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Further, from the Bethe ansatz result Ijllfl we have ph 
1/2-k. Thus the distribution function is obtained as 



n(k) = 2np(k) 



1 



e (k2-p,)/K B T _ l 



(15) 



which represents the Bose statistics with a = 0. This 
relation is easily seen from the GES result (|13|l with a = 
0. 

In the strong coupling regime and at low tempera- 
tures, the relation between the densities p and ph is 
2ir (ap + ph) ~ 1 with nonmutal GES determined by 
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The density relation gives the most probable distribution 
n(e) = 2np to be [Il| 



n(e) 



w(e) 



a In (1 — w(e)) — lnw(e) 
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(17) 
(18) 



The relations H17fl and (|18fl provide an alternative path 
to the thermodynamics of the ID interacting Bose gas 
via well established results in quan tum statistics. In par- 
ticular, following Isakov et al. ]E3 and Iguchi [&|, at low 
temperatures, i.e., for T < T<j, where Td — is the 

quantum degeneracy temperature, the thermodynamics 
of ideal particles obeying GES can be derived from the 
fractional statistics distribution n(e) via Sommerfeld ex- 
pansion. Thermodynamic properties such as the total 
energy E and the pressure may be calculated directly 
from 1)20(1 in terms of the integral I[f] given in Ref. 
for ideal particles obeying GES. 

In order to calculate the thermodynamic properties, 
we write the particle number and total energy as 



N = / deG(e)n(e), (19) 
Jo 

/>oo 

E = / deG{e)n{e)e, (20) 
in terms of the density of states G(e) given by 



G(e)=L/[2 nl -e 



r, 2 



(21) 



In this way we find the chemical potential |4lj 

p = p [1 + c 2 t 2 + c 3 t 3 + c 4 i 4 + 0(t 5 )] , (22) 

in terms of the effective temperature t — KbT/pq, where 
po = ■^K 2 n 2 a 2 . The first few coefficients are given by 



C2 
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c 3 = --C(3)a(l - a), 



c 4 = — — -a (3 - 2a + 3a 2 ) , 
144 v / 



(23) 
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with the constant C(3) = J2n=i V™ 3 - 1-20206. 

On the other hand, for the strongly interacting Bosc 
gas, the distribution function obtained from the TBA 
results ((Tr) and ljT2l is 



n(e) = 2itp(t) 



a(l 



(24) 



with the dispersion relation e = i—k 2 . Here the chemical 

1 Am 

potential is given by 



A* = A*o 



1 



(25) 



It follows that in the strong coupling limit the chemical 
potential ill' 21) derived from GES coincides with the TBA 
result 125(1 . For 7 ^> 1 they coincide to the order shown, 
with the coefficient C3 in (GES-mu) vanishing. For com- 
parison, the distribution profiles obtained from GES l|17|) 
with (UBJ and the TBA result JUJ with (J55J) are shown 
in Figure ^ There is good agreement between the re- 
sults obtained from the two approaches for strongly in- 
teracting bosons. At low temperatures the distributions 
for strongly interacting bosons clearly deviate from pure 
Fermi statistics. As seen in Figure ^ the height of the 
distribution function n(e) approaches the Fermi statis- 
tics value n(e) = 1 with increasing interaction strength 
7 at zero temperature. During the fermionization pro- 
cess, i.e. as 7 increases from large but finite values to 
infinity, more than one particle is allowed to occupy a 
single quantum state. This signature is attributable to 
the collective behaviour of interacting bosons. 



At low temperatures, the thermodynamic properties 
derived from GES coincide with the ones derived from 
the TBA H3- The total energy obtained from the TBA 
result J23J is 
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in terms of the effective temperature r = KbT/T^ with 



degeneracy temperature Td 



2m ' 



and where En = 



^N/jq is the ground state energy © at zero temperature. 
Figure |21 shows a plot of the total energy per particle in 
units of the Fermi energy 



E F = 



2m~3' 



2 2 
in, 



(27) 



as a function of the coupling constant 7 and the effective 
temperature r. It is apparent that the energy increases 
slowly as the effective temperature r increases, with the 
total energy mainly dominated by the dynamical inter- 
action. The thermal fluctuations are largely suppressed. 
The energy curve approaches that obtained with Fermi 
statistics as 7 — > 00, namely in the TG limit. 
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FIG. 1: Comparison between the GES and TBA results for 
the most probable distribution profiles n(e) for the values 
7 = 10 (left panel) and 7 = 20 (right panel) at different values 
of the degeneracy temperature r = KbT /Td- At zero tem- 
perature n(e) = 1/q which leads to a Fermi surface at e = fin- 
This surface gradually decreases with increasing temperature 
due to a large number of holes appearing below the surface. 
Pure Fermi statistics are obtained in the limit 7 — » 00. The 
solid lines show the most probable GES distribution for ideal 
particles derived from 1171 with 1221 . The symbols indicate 
the corresponding distributions evaluated from the TBA re- 
sult 112 11 for interacting bosons. 



The specific heat is given by = (fip)^ L with result 
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Figure [3] shows a plot of the specific heat as a func- 
tion of the effective temperature for different interaction 
strength. For small temperatures the specific heat is al- 
most linearly increasing with r. Similar behaviour was 
found in the (/-deformed Bose gas jH^]. For strongly in- 
teracting bosons the specific heat deviates from the free 
Fermi curve as 7 decreases. This is mainly because the 
dynamical interaction 7 lowers the entropy in fermion- 
ization. In general the specific heat reveals an important 
signature of the quantum statistics of interacting many- 
body systems as it sensitively depends on the GES pa- 
rameter a. 



We now consider the local pair correlations for the 
strongly interacting ID Bose gas and their role in re- 
vealing GES. The one-body and two-body correlation 
functions are key quantities for describing the quan- 
tum de gen erate signature of ID interacting systems 
|29l l30ll3ll I32I |U The momentum distribution and 
the static structure factor can be obtained from the 
Fourier transform of the one-body and two-body corre- 
lation functions, respectively. The local three-body cor- 
relation gives the three-body recombination rate. Uni- 
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FIG. 2: The total energy per particle in units of the free Fermi 
energy _Ef as a function of the effective temperature r and the 
coupling strength 7. 
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FIG. 3: Specific heat (H§ in units of NK B as a function of 
the degeneracy temperature r for the interaction parameter 
values 7 = 20, 30, 50, 100, 200, 00. The values of the statistical 
parameter a vary with respect to 7 according to 1161 . 



versality in correlation functions is of particular inter- 
est in theory HI HI El El El El and experiment 
0,13- In particular, the local pair correlations for the ID 
Bose gas have been observed experimentally by measur- 
ing photoassociation rates [8|. Physically, the local pair 
correlation is a measure of the probability of observing 
two particles in the same place. 

In the grand canonical description, the two-particle lo- 
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FIG. 4: The local pair correlation <? (0) as a function of the 
effective temperature r and the coupling strength 7. 



cal correlation is given by [T2I l2fj l8ll | 



(29) 



where f = fi— 2E/N is the free energy per particle. In 
this way the result 



fJ 
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for the two-particle local correlation at low temperatures 
follows from the TBA. Figure 0] shows the local pairing 
correlation as a function of the effective temperature r 
and the interaction strength 7. We see that the dynam- 
ical interaction dramatically reduces the pair correlation 
due to decoherence between individual wave functions. 
It tends to zero as 7 — > 00, namely in the TG limit. 
Moreover, the pair correlation is only weakly dependent 
on the temperature because the thermal fluctuations are 
suppressed at the temperatures T <?C Th. The experimen- 
tal values for the local pair correlations y| are in excellent 
agreement with the theoretical result H30JI for over a wide 
range of interaction strength for r f» 0. 



III. DISCUSSION 

We have demonstrated that the statistical profiles of 
the strongly interacting ID Bose gas at low temperatures 
are equivalent to those of a gas of ideal particles obeying 
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nonmutual GES. For the strong coupling limit 7 — » 00, 
the TG gas, the equivalence is with Fermi statistics. The 
deviation from Fermi statistics during the fermionization 
process is described by the relation a ~ 1 — 2/7, where a 
is a measure of the GES. In the strong coupling regime 
the local pair correlation function (|30(l may also be writ- 
ten as 

This suggests the possibility of observing the quantum 
degenerate behaviour of ideal particles obeying GES 
through experiments on ID interacting systems. In this 
sense the recent experimental measurement of the local 
pair correlation function in a ID Bose gas of 87 Rb atoms 
also provides a measurement of GES. In particular, 



the deviation of the GES parameter a away from the 
pure Fermi statistics value a = 1. This gives a quanti- 
tative description of the fermionization process. Other 
thermodynamic properties, such as the distribution pro- 
files depicted in Figure 1 and the specific heat, are also 
sensitive to the statistics. They also thus provide a way of 
exploring fractional statistics in the strongly interacting 
ID Bose gas. 
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